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Abstract 

We calculate the 2nd and 4th order quark number susceptibihties at zero baryochemical potential, using a PNJL approach 

- - -and an approach which includes, in a single model, quark and hadronic degrees of freedom. We observe that the 

. susceptibilities are very sensitive to possible quark-quark vector interactions. Compared to lattice data our results 

f~^ 'suggest that above Tc any mean field type of repulsive vector interaction can be excluded from model calculations. 

CNJ Below Tc our results show only very weak sensitivity on the strength of the quark and hadronic vector interaction. The 

^ best description of lattice data around Tc is obtained for a case of coexistence of hadronic and quark degrees of freedom 
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Many recent experimental programs in heavy ion physics, 
'e.g. at the relativistic heavy ion collider and the planned 
FAIR facility, are aimed at a better understanding of bulk 
properties of QCD matter. In particular the experimental 
confirmation of the deconfinement and chiral phase transi- 
tions as well as the search for the critical end point (CeP) 
'are of great interest to the community. Lattice results 
■at ^B = suggest that both phase transitions exhibit a 
smooth crossover. Suffering from the so called sign prob- 
lem at finite baryo chemical potential, lattice studies have 
so far not been able to constitute a consistent picture of the 
■phase diagram and the existence of the CeP [l|, |2| ■ While 
a Taylor expansion of lattice results at /if, = predict the 
existence of a critical end point these results may strongly 
depend on the order to which the Taylor coefHcients can be 
Evaluated [3|. Lattice calculations at imaginary chemical 
potential on the other hand show no indication for a CeP, 
though these calculations are done on on a coarse lattice 
,with a fermion action that is known to have large dis- 
cretization errors. Other studies trying to locate the CeP 
often rely on the applicability of various effective models 
where CeP's existence and location seems very sensitive 
on the vector coupling strength at and below Tc [J]. 

1. The PNJL model 

The PNJL model was introduced in [sl |6| as an ef- 
fective chiral quasi-quark model that incorporates a mean 
field like coupling to a color background field. It has often 
been shown to reproduce many general features of lattice 
results at /is = 0, i, i, [13, [ill [H, [H, Q E El, [ll, 
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i, [H, [H, [13, [m, [H, [H, [13, [H, [H, [13, [2^ . In our com- 
parison we will use a very basic parametrization of the 
two-flavor PNJL model and extend it to incorporate a re- 
pulsive vector interaction. The thermodynamic potential 
of our parametrization reads: 



n = C/($, $*, T) + ctVSGs - a;V2Gv 
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Here the grand canonical potential includes contributions 

of sates with 1, 2 and 3 times the single quark energy. Note 

that the 3 quark state does not couple to the Polyakov 

loop. 

The (traced) Polyakov loop <I> was introduced as: 



$ = 1/3 Tr e"^/^ 



(3) 



with (j> = A4, a. background color gauge field. The ther- 
modynamics of $ (and <&*) are controlled by the effective 
potential t/($,$*,T) Q: 



U = --a(T)$$* 

+ b{T) ln[l - 6$$* -h 4($3$*3) - 3($$*)2] (4) 



Preprint submitted to Elsevier 



January 6, 2011 



with a{T) = aoT^ + aiT^T^ + a2T^T^ 



b{T) = hT^T. 



This choice of effective potential satisfies the Z{3) cen- 
ter symmetry of the pure gauge Lagrangian. In the con- 
fined phase, U has a minimum at $ = 0, while above the 
critical Temperature Tq its minimum is shifted to finite 
values of $. The logarithmic term appears from the Haar 
measure of the group integration with respect to the SU(3) 
Polyakov loop matrix. The parameters ag, oi, a2 and 63 are 
fixed, as in [29|, by demanding a first order phase transi- 
tion in the pure gauge sector at Tq = 270MeV, and that 
the Stefan-Boltzmann limit is reached for T — > 00. 

The dynamical mass of the quarks m — niQ ~ a — 
"^0 ^ Gs (^^) is the same as in the NJL model and the 
vector coupling induces an effective chemical potential for 
the quarks fi* = fiq + uj = Hq + Gy (^t*). The two aux- 
iliary fields a and w are controlled by the potential terms 
and the last term includes the difference AEp between the 
quasi particle energy and the energy of free quarks. The 
NJL part of the model has 4 parameters, the bare quark 
mass for the u- and d-quarks (assuming isospin symmetry), 
the three-momentum cutoff of the quark-loop integration 
A and the coupling strengths Gs and Gy- To reproduce 
realistic values for the pion mass and decay constant as 
well as the chiral condensate, we take these values to be 
@]: mu,d = 5.5 MeV, Gs = 10.08 GeV"^A = 651 MeV 
{Gv will be left as a model parameter to study the influ- 
ence of the vector coupling on our results) . 
The self consistent solutions are obtained by minimizing 
the thermodynamic potential with respect to the fields a, 
uj, $ and $*. 

2. The QH model 

To estimate the infiuence, of hadronic contributions, 
on the susceptibilities we will compare the results from 
the PNJL model with those obtained from a model where 
quark and hadronic degrees of freedom are present in a 
single partition function (Quark-Hadron model). 
In the following we will shortly describe the different com- 
ponents of the model, for a more detailed discussion we 
refer to [SO*. The hadronic part is described by a flavor- 
SU(3) model, which is an extension of a non-linear repre- 
sentation of a cr-uj model including the pseudo-scalar and 
vector nonets of mesons and the baryonic octet [3l|, |32|, |33[ . 
Besides the kinetic energy term for hadrons and quarks, 
the terms: 

Ltnt = - Li i^ihoigiL^^ + 9t<li4') + "^*]V'^, (5) 
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represent the interactions between baryons and vector and 
scalar mesons, the self-interactions of scalar and vector 
mesons, and an explicitly chiral symmetry breaking term. 
The index i denotes the baryon octet and the different 
quark flavors. Here, the mesonic condensates (determined 
in mean-field approximation) included are the vector-isoscalars 
UJ and (/), and the scalar-isoscalars a and C, (strange quark- 
antiquark state) . The last four terms of eqn. © were intro- 
duced to model the QCD trace anomaly, where the dilaton 
field X can be identified with the gluon condensate. 

The effective masses of the baryons and quarks are gen- 
erated by the scalar mesons except for an explicit mass 
term {Sttin = 120 MeV, Sruq = 5 MeV and Sm^ = 105 
MeV for the strange quark), m* — gia-cr + gi(^( + Srui. 
Vector type interactions introduce an effective chemical 
potential for the quarks and baryons, generated by the 
coupling to the vector mesons: fj,* — ^i — ft^w — gi^t/i. 
The coupling constants for the baryons [3J] are chosen to 
reproduce the vacuum masses of the baryons, nuclear sat- 
uration properties and asymmetry energy as well as the A- 
hyperon optical potential. The vacuum expectation values 
of the scalar mesons are constrained by reproducing the 
pion and kaon decay constants. For the quarks we chose 
the following coupling parameters: gq„ — gsQ — 4.0, while 
the vector coupling strength gq^^ is left as free parameter. 
The coupling of the quarks to the Polyakov loop is intro- 
duced in the thermal energy of the quarks 30] . All ther- 
modynamical quantities, energy density e, entropy density 
s as well as the densities of the different particle species 
Pi , can be derived from the grand canonical potential: 



' ' — i-iint ^ J^ri 



^th - u 
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Here Clth includes the heat bath of hadronic and quark 
quasi particles. The effective potential [/(<!>, $*,T) which 
controls the dynamics of the Polyakov-loop has the form 
of eqn.dl]). 

Since we expect the hadronic contribution to disappear, at 
least at some point above Tc , we included effects of finite- 
volume particles to effectively suppresses the hadronic de- 
grees of freedom, when deconfinement is achieved. Includ- 
ing these effects in a thermodynamic model for hadronic 
matter, was proposed some time ago [33, 'Sa, ISll |38| ■ We 
will use an ansatz similar to that in [3 ft . i4Q] , but modify it 
to also treat the point like quark degrees of freedom con- 
sistently. We introduce the chemical potential Jli which 
is connected to the real chemical potential /i*, of the i- 
th particle species by the relation: fli ~ jjl* ~ Vi P {P is 
the sum over all partial pressures and Vi the volume of a 
hadron). All thermodynamic quantities can then be calcu- 
lated with respect to the temperature T and the new chem- 
ical potentials Jli. To be thermodynamically consistent, all 
densities have to be multiplied by a volume correction fac- 
tor /, which is the ratio of the total volume V and the 
reduced volume V , not being occupied. In this configura- 
tion the chemical potentials of the hadrons are decreased 
by the quarks, but not vice versa. As the quarks start ap- 
pearing they effectively suppress the hadrons by changing 
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Figure 1: The second order quark number susceptibility from the 
PNJL model, with different strengths of the vector interaction, as a 
function of T over Tc- Black solid line: no vector interaction, red 
dashed line: Gy = Gs/2, green dotted line: Gy = Gg, blue dash 
dotted Une Gy = "^Gs- 



their chemical potential, while the quarks are only affected 
through the volume correction factor /. Our description 
of the excluded volume effects is admittedly simplified and 
parameter dependent, but it enables us to describe a phase 
transition from hadronic to quark degrees of freedom in a 
quite natural and thermodynamically consistent manner. 

As has been mentioned above, the Lagrangian of the 
chiral model contains dilaton terms to model the scale 
anomaly. These terms constrain the chiral condensate, if 
the dilaton is frozen at its ground state value xo- On the 
other hand, as deconfinement is realized, the expectation 
value of the chiral condensate should vanish at some point. 
On account of this we will couple the Polyakov loop to the 
dilaton in the following way: 



X = xo (1 - 0.5(<i><i>*)) 



(8) 



Assuming a hard part for the dilaton field which essentially 
stays unchanged and a soft part, which vanishes when de- 
confinement is realized. Hence, allowing the chiral conden- 
sate to also approach zero. 

For a more detailed discussion of the model and compar- 
isons with lattice data we refer to [30[. 

3. Susceptibilities 

Lattice results at finite chemical potentials are often 
obtained as Taylor expansion of the thermodynamic quan- 
tities in the parameter ^/T around zero chemical potential 
4l|. In the Taylor expansion of the pressure p = — fi, the 
coefficients, which can be identified with the quark number 
susceptibilities, follow from: 
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Figure 2: The second order quark number susceptibility from the 
QHM model, with different strengths of the vector interaction, as a 
function of T over T^. Black solid line: gjvij 7^ Oi 9qt^ = 0, i; = 1, red 
dashed line: g^i^j 7^ 0, gq^j = gNi^/3, v = 1, green dotted line: gNu = 
0, gquj = 0, V = 1, blue dashed dotted line: gjvij 7^ 0, gq^j = 0, i) = 0, 
orange dash dot dot line: gjvuj = 0, gqui = 0, ti = 0. 



Cn{T) = 



i5"(p(r,MB)/T4) 



n\ difiB/ry 



(10) 



fJ.B=0 
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In our approach we explicitly calculate the pressure at 
finite /is and then extract the expansion coefficients nu- 
merically. The results for the second coefficient calculated 
for the PNJL model and QH model, compared to lattice 
results [43, are shown in figures ^ and ^. One can 
clearly observe that the best description can be obtained 
when quark vector interactions are turned off. Any form 
of repulsive interaction strongly decreases the value of the 
second order coefficient above Tc- The lattice results on 
the other hand quickly reach a value that is expected for a 
non-interacting gas of quarks, even right above Tc- This is 
in fact surprising as other thermodynamic quantities tend 
to favor a picture with a wide region around Tc where in- 
teractions are strong. This behavior was also reproduced 
by the PNJL model and the QH model while both fail to 
describe the fast increase in C2 right above Tc- 

Next we try to disentangle the hadronic contribution 
to the second coefficient. As can be seen in figure (O 
the dependence on the strength of the repulsive interac- 
tion below Tc is rather small and one can not exclude any 
scenario. In the crossover region, around T^ differences 
become obvious. 

The solid black line in Figure (I2|) displays the result for 
C2 using the standard parametrization of the QH model as 
described in [30| without any repulsive quark-quark inter- 
actions. The value of C2 only slowly approaches 1, which 
is mainly due to the fact that the value of the chiral con- 
densate drops to rather slow in our model and therefore 
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Figure 3: The forth order quark number susceptibihty from the 
PNJL model, with different strengths of the vector interaction, as 
a function of T over Tc- Black solid line: no vector interaction, red 
dashed line: Gy = Gs/2, green dotted line: Gy = Gg, blue dash 
dotted hne: Gv = 2Gs- 



the quark masses do not decrease as fast as in the PNJL 
approach. 

Including repulsive vector interactions for the quarks 
gives a result which is similar to the one obtained from 
the PNJL (red dashed line). Here the repulsive interac- 
tions strongly decrease the value of C2 . 

Turning off all repulsive interactions, vector interac- 
tions for quarks and hadrons as well as the excluded vol- 
ume corrections (orange dash dot dotted line), leads to 
an drastic overestimation of C2. This is expected as all 
hadronic degrees of freedom are present at and above Tc if 
the excluded volume effects are turned off. Therefore one 
largely overestimates the effective degrees of freedom. 

On the other hand if the repulsive vector interactions 
are turned on only for the hadrons one obtains a rather 
good description of the lattice results (blue dash dotted 
line). In this parametrization the hadrons are also still 
present in the system up to arbitrary high temperatures, 
as the excluded volume effects are still turned of, and 
therefore all thermodynamic quantities are largely over 
predicted. 

To remove the hadronic contributions from the system 
we introduced excluded volume corrections as described 
above. This leads to a pronounced dip in C2 above Tc 
(green dotted line), indicating that our excluded volume 
approach is either to simplified or all hadronic contribu- 
tions are already vanishing completely right above Tc or 
that the lattice results are still not accurate enough to suf- 
ficiently resolve effects of hadron hadron interactions. 

The fourth order coefficients calculated for the PNJL 




Figure 4: The fourth order quark number susceptibility from the 
QHM model, with different strengths of the vector interaction, as a 
function of T over Tc- Black solid line: gjvu ¥" Oi Sq" = 0, i" = 1, red 
dashed line: gjvt^ ^ 0, gq^ = gNi^/3, v = 1, green dotted line: gjvt^ = 
0,gquj = 0,^' = 1, blue dashed dotted line: gjvuj ¥" ^^9qt^ = 0, i) = 0, 
orange dash dot dot line: gisiui = 0, gqu, = 0,v = 0. 



and QH model are shown in figures [3] and SI Although the 
errors on the lattice results are still significant our results 
for C4 support the statements made for C2- At tempera- 
tures above Tc a gas of free non interacting quarks, with- 
out any hadronic contribution, gives the only reasonable 
description of the data already slightly above Tc- While 
the height of the peak does depend on the strength of the 
vector coupling, we observe also a strong dependence on 
the slope of the change of the order parameters around Tc, 
which is larger in the QH model as compared to the PNJL 
model. Notice that the parametrization shown as the blue 
dash-dotted line, which gave the best description for C2, 
overestimates C4 around Tc drastically. This is simply be- 
cause all thermodynamic quantities are over predicted in 
this case and therefore the densities as well as the order 
parameters increase steeper as they do in the lattice cal- 
culations. 

4. Conclusion 

We have compared results for the quark number sus- 
ceptibilities at /is = obtained from two mean field mod- 
els to recent results from lattice calculations. Both models 
strongly indicate that the EoS of QCD above Tc seems to 
be composed purely of a gas of non-interacting quasi par- 
ticles. This finding is in agreement with recent work by 
[43[ where the extracted quark vector coupling strength 
approaches zero very fast around Tc- On the other hand, 
lattice observables like the interaction measure and the 
normalized Polyakov loop indicate a large region above 
the critical temperature where the hot QCD medium is far 



from being an ideal gas. Around T^ repulsive hadronic in- 
teractions are supported by our results with the QH model. 
To describe the steep increase of C2 around Tc we need to 
introduce hadronic contributions up to right above T^. At 
even larger temperature the hadrons seem to be replaced 
by a almost non interacting gas of quarks. This offers an 
intriguing implication concerning the CeP. In mean field 
studies with the PNJL [4| as well as in the HQ model [30| , 
a general feature of including repulsive interactions is that 
the CeP is moved to larger chemical potentials or even 
disappears completely. From this point of view, a vanish- 
ing quark vector interaction would favor the existence of 
a CeP. 

On the other hand, the phase structure of QCD up to Tc 
and likely even slightly above it could be determined by 
hadronic interactions. In model calculations without any 
repulsive quark interactions the quarks often appear at 
rather small chemical potentials (at T = 0) making a rea- 
sonable description of the nuclear ground state difficult. 
A repulsive hadronic vector interaction, which is required 
on order to reproduce the properties of a saturated nuclear 
ground state, may move the CeP to larger chemical poten- 
tials or even remove it completely from the phase diagram 
in accordance to mean- field results [J, |43j . 
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